NUWC-NPT  Technical  Report  11,762 
1  September  2006 


Bending  Behavior  of  Plain-Woven  Fabric 
Air  Beams:  Fluid-Structure  Interaction 
Approach 

Paul  V.  Cavallaro 

Ranges,  Engineering,  and  Analysis  Department 

Ali  M.  Sadegh 

The  City  College  of  New  York 

Claudia  J.  Quigley 

U.S.  Army  Research,  Development  and  Engineering  Command 


W&SEA 

WARFARE  CENTERS 


NEWPORT 

Naval  Undersea  Warfare  Center  Division 

Newport,  Rhode  Island 


Approved  for  public  release;  distribution  is  unlimited. 


PREFACE 


This  report  was  funded  under  NUWC  Division  Newport 
Assignment  Number  TD0207,  principal  investigator  Paul  V. 
Cavallaro  (Code  70T).  The  research  was  conducted  in  support  of 
the  Center  of  Excellence  for  Inflatable  Composite  Structures  at  the 
U.S.  Army  Natick  Soldier  Center  (NSC)  in  Natick,  MA,  through 
Military  Interagency  Purchase  Request  Number  6CS6R003 19. 

The  NSC  project  officer  is  Jean  Hampel. 

The  technical  reviewer  for  this  report  was  Andrew  J.  Hull 
(Code  8212). 

The  authors  especially  thanJk  Karen  Horak  of  the  NSC  for 
coordinating  technical  discussions  and  information  exchanges  with 
industry,  academia,  and  other  DoD  agencies;  Joseph  J.  Stace  (Code 
1512)  for  conducting  the  multi-axial  fabric  tests  and  extensive  data 
reductions;  and  Adam  Turner  of  the  University  of  Maine  at  Orono 
for  providing  the  high  quality  fabric  material  samples  that  were 
crucial  to  obtaining  exceptional  experimental  results. 


Reviewed  and  Approved:  1  September  2006 


Harriet  L.  Coleman 

Head,  Ranges,  Engineering,  and  Analysis  Department 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMB  No.  0704-0188 


Public  reporting  for  this  collection  of  information  Is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  Instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  Information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of 
information,  including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway, 
Suite  1204,  Arlington,  VA  22202-4302,  and  to  the  Office  of  Management  and  Budget,  Paperwork  Reduction  Project  (0704-0188),  Washington,  DC  20503. 


1 .  AGENCY  USE  ONLY  (Leave  blank)  2.  REPORT  DATE  3.  REPORT  TYPE  AND  DATES  COVERED 

1  September  2006 


4.  TITLE  AND  SUBTITLE 

Bending  Behavior  of  Plain-Woven  Fabric  Air  Beams:  Fluid-Structure  Interaction 
Approach 


6.  AUTHOR(S) 


5.  FUNDING  NUMBERS 


Paul  V.  Cavallaro,  Ali  M.  Sadegh,  Claudia  J.  Quigley 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Naval  Undersea  Warfare  Center  Division 
1176  Howell  Street 
Newport,  Rl  02841-1708 


9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

U.S.  Army  Natick  Soldier  Center 

Kansas  Street 

Natick,  MA  01760-5018 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 

TR  11,762 


10.  SPONSORING/MONITORING 
AGENCY  REPORT  NUMBER 


12a.  DISTRIBUTION/AVAILABILITY  STATEMENT  12b.  DISTRIBUTION  CODE 

Approved  for  public  release;  distribution  is  unlimited. 


13.  ABSTRACT  (Maximum  200  words) 

A  swatch  of  plain-woven  fabric  was  subjected  to  biaxial  tests  to  determine  its  material  characteristics.  The  stress-strain 
relationships  of  the  fabric  were  determined  and  used  directly  in  finite  element  models  of  an  air  beam  that  was  assumed  to  be 
constructed  of  the  same  fabric  and  subjected  to  inflation  and  bending  events.  The  structural  responses  to  these  events  were 
obtained  using  the  ABAQUS/Explicit  Finite  element  solver  for  a  range  of  pressures,  including  those  considered  to  be  typical  in  safe 
operations  of  air-inflated  structures.  The  models  accounted  for  the  flu  id -structure  interactions  between  the  air  and  the  fabric.  The 
air  was  treated  as  a  compressible  fluid  in  accordance  with  the  Ideal  Gas  Law  and  was  subjected  to  adiabatic  constraints  during 
bending.  The  fabric  was  represented  with  membrane  elements.  Several  constitutive  cases,  including  linear  elasticity  and 
hyperelasticity,  were  studied.  The  bending  behavior  for  each  constitutive  case  is  presented,  followed  by  a  discussion  of  the  use  and 
limitations  of  the  cases. 


14.  SUBJECT  TERMS 

Inflated  Structures 
Finite  Element  Analysis 


17.  SECURITY  CLASSIFICATION 
OF  REPORT 

Unclassified 


Air  Beams 
Plain-Woven  Fabric 


18.  SECURITY  CLASSIFICATION 
OF  THIS  PAGE 

Unclassified 


Crimp  Interchange 
Explicit  Methods 


19.  SECURITY  CLASSIFICATION 
OF  ABSTRACT 

Unclassified 


15.  NUMBER  OF  PAGES 

34 


16.  PRICE  CODE 


20.  LIMITATION  OF  ABSTRACT 


NSN  7540-01-280-5500 


Standard  Form  298  (Rev.  2-89) 
Prescribed  by  ANSI  Std.  Z39-18 
298-102 


TABLE  OF  CONTENTS 


Page 

LIST  OF  TABLES . ii 

LIST  OF  ABBREVIATIONS  AND  ACRONYMS . iii 

INTRODUCTION . 1 

DESCRIPTION  OF  PLAIN-WOVEN  FABRIC . 1 

FABRIC  STRESS-STRAIN  BEHAVIOR— EXPERIMENT . 3 

AIR  BEAM  FINITE  ELEMENT  MODELS . 8 

ANALYSIS  AND  RESULTS . 10 

Case  1 — Linear  Elasticity,  Constant  Modulus  for  All  Pressures . 10 

Case  2 — Linear  Elasticity,  Use  of  Instantaneous  Tangent  Moduli . 12 

Case  3 — Hyperelasticity . 14 

DISCUSSION . 17 

CONCLUSIONS . 18 

REFERENCES . 19 

APPENDIX  A  —  IDEAL  GAS  EQUATION  OF  STATE  (EOS) . A-l 

APPENDIX  B  —  REDUCED  POLYNOMIAL  (N=  3)  STRAIN  ENERGY  POTENTIAL . B-l 

LIST  OF  ILLUSTRATIONS 

Figure  Page 

1  Image  of  Plain-Woven  Polyester  Fabric . 2 

2  Stress  State  in  Cylinders  Subjected  to  Internal  Pressure . 3 

3  Combined  Multi-Axial  and  Shear  Test  Fixture  (U.S.  Patent  No.  6,860,156) . 3 

4  Dimensions  of  Fabric  Swatch  Specimen . 4 

5  Plain-Woven  Polyester  Fabric  Specimen  Subjected  to  Multi-Axial  and  Shear 

Tests . 4 


l 


LIST  OF  ILLUSTRATIONS  (Cont'd) 


Figure  Page 

6  Uniaxial  and  Biaxial  Stress-Strain  Curves  for  Plain-Woven  Polyester  Fabric . 5 

7  Stress  Ratio  S  vs  Equi-Biaxial  Strain  e . 5 

8  Etan  vs  crWarP  Measured  in  the  Equi-Biaxial  Extension  Test . 6 

9  r  vs  y  Plotted  with  P  at  10  psi  ( awarp  =  637  psi) . 6 

10  G  vs  y  with  P  at  10  psi  (owarp  =  637  psi) . 7 

1 1  Kinematics  of  Shear  Deformations  in  Uncoated  Plain-Woven  Fabrics . 7 

1 2  Air  Beam  Finite  Element  Model  Subjected  to  Four-Point  Bending . 8 

1 3  Transverse  Shear  Deformations  During  Bending . 9 

14  fVex,,  PdV-work,  and  A Es,ram  vs  P  for  a  Linearly  Elastic  Fabric — Case  1 . 1 1 

1 5  F 'reaa  vs  dmid  Curves  for  Linear  Elastic  Material  Model  with  P  at  1 0  to  60  psi — 

Case  1  . 11 

16  Wext,  A Esiram,  and  PdV-work  vs  P  for  a  Linearly  Elastic  Fabric  When  E',an  Used 

with  P  at  10,  15,  and  20  psi — Case  2 . 13 

17  F react  vs  bmut  for  a  Linearly  Elastic  Fabric  When  E',an  Used  with  P  at  10,  15,  and 

20  psi — Case  2 . 13 

1 8  Wext,  strain i  and  PdV-work  vs  P  for  a  Hyperelastic  Fabric:  Reduced  Polynomial 

(N=  3)  Strain  Energy  Potential  Applied  with  -Pat  10,  15,  and  20  psi — Case  3 . 15 

1 9  F reaa  vs  Smjd  for  Hyperelastic  Fabric  Model  with  P  at  10,  15,  and  20  psi — Case  3 . 15 

20  Time  History  Plot  of  owarp  for  Inflation  and  Bending  Steps  with  P  at  10,  15,  and 

20  psi — Case  3 . 16 

2 1  Plots  of  owarp  and  owef,  vs  Strain  for  P  at  10,  15,  and  20  psi — Case  3 . 16 


LIST  OF  TABLES 


Table  Page 

1  Details  of  the  Plain-Woven  Polyester  at  the  Yam  and  Fabric  Levels . 2 

2  Summary  of  Pressure,  Volume,  and  Energy  Changes  During  Bending  Step — 

Case  1 . 12 

3  Owarp  and  E',an  as  Functions  of  P . 12 

4  Summary  of  Pressure,  Volume,  and  Energy  Changes  During  Bending  Step — 

Case  2 . 14 

5  Summary  of  Pressure,  Volume,  and  Energy  Changes  During  Bending  Step  for  the 

Hyperelastic  Fabric  Model — Case  3 . 15 


ii 


LIST  OF  ABBREVIATIONS,  ACRONYMS,  AND  SYMBOLS 


Si  warp 
A  weft 

c 

Co 

A 

E 

F' 

tan 
Ed  is 

Ejnt 

E kinetic 

EOS 

Estrain 
Eton 
E react 

G 

l 

f 

Efabric 

Eyarn 

hdw 

Ef wrinkle 

N 

NSC 

NUWC 

P 

Pa 

PdV-work 

r 

R 

R 

5 
S° 
u 

V 

wex, 

YDR 

Y 
A 

&load _pt 
&mid 

e 

<r 

6 


Warp  yam  cross-sectional  area 
Weft  yam  cross-sectional  area 
Crimp  content  (%) 

Temperature  dependent  material  parameter 

Temperature  dependent  material  parameter 

Elastic  modulus 

Instantaneous  tangent  modulus 

Viscous  dissipation  energy 

Internal  energy 

Kinetic  energy 

Equation  of  state 

Strain  energy 

Tangent  modulus 

Sum  of  load  point  reaction  forces 

Shear  modulus 

First  deviatoric  strain 

Elastic  volume  ratio 

Reference  length  of  yam  in  fabric 

Length  of  yam  extracted  from  fabric 

Molecular  weight 

Wrinkling  moment 

Material  parameter 

U.S.  Army  Natick  Soldier  Center 

Naval  Undersea  Warfare  Center 

Inflation  pressure 

Ambient  pressure 

Air  compressibility 

Radius  of  cylinder 

Gas  constant 

Universal  gas  constant 

Stress  ratio 

Converged  stress  ratio 
Strain  energy  per  volume 
Air  volume 
External  work  done 
Yam  density  ratio 
Shear  strain 

Change  between  bending  and  inflation  end  states 

Load  point  displacement 

Mid-span  deflection 

Strain  due  to  equi-biaxial  extension 

Crimp  ratio 

Current  temperature 


LIST  OF  ABBREVIA  TIONS,  ACRONYMS,  AND  SYMBOLS  (Cont’d) 


9Z 

Absolute  zero  temperature 

Mo 

Initial  hyperelastic  shear  modulus 

V 

Poisson’s  ratio 

P 

Density  of  air 

&warp 

Warp  direction  stress 

Oweft 

Weft  direction  stress 

T 

Shear  stress 

IV 


BENDING  BEHAVIOR  OF  PLAIN-WOVEN  FABRIC  AIR  BEAMS: 
FLUID-STRUCTURE  INTERACTION  APPROACH 

INTRODUCTION 
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This  report  extends  previous  experimental  and  numerical  modeling  efforts  ’  that  addressed 
the  bending  behavior  of  uncoated,  plain-woven  fabric  air  beams.  The  prior  efforts  showed  that 
the  elastic  and  shear  moduli  were  independent  of  the  elastic  modulus  of  the  yams  but  were 
dependent  upon  inflation  pressure,  fabric  construction,  and  crimp  interchange.2  The  models  also 
treated  the  stiffening  influence  of  inflation  pressure  on  only  the  constitutive  properties  of  the 
fabric.  Changes  in  pressure  and  volume  were  not  considered. 

The  objectives  of  the  current  research  were  to  (1)  characterize  the  mechanical  properties  of  a 
structural  fabric  through  swatch-level  tests  by  using  a  combined  biaxial  tension  and  shear  test 
fixture,3  (2)  input  the  resulting  biaxial  stress-strain  behavior  directly  into  air  beam  finite  element 
models,  (3)  predict  the  bending  behavior  of  the  air  beam  while  accounting  for  fluid-structure 
interactions  (i.e.,  air  compressibility),  (4)  identify  the  limitations  of  using  linear  elasticity  and 
hyperelasticity  to  represent  the  fabric  constitutive  behavior,  and  (5)  assess  the  contributions  of  air 
compressibility  (as  a  form  of  PdV- work)  and  fabric  strain  energy  on  air  beam  bending  behavior. 

A  swatch  of  plain-woven  fabric  was  subjected  to  biaxial  tests  to  characterize  its  mechanical 
properties.  The  stress-strain  relations  of  the  fabric  were  determined  and  used  directly  in  finite 
element  models  of  an  air  beam  that  was  assumed  to  be  constructed  with  the  same  fabric.  The 
structural  responses  to  inflation  and  bending  events  were  obtained  using  the  ABAQUS/Explicit4 
finite  element  solver  for  a  range  of  pressures,  including  those  considered  to  be  typical  in  safe 
operations  of  air-inflated  structures.  The  models  accounted  for  the  fluid-structure  interactions 
between  the  air  and  the  fabric.  The  air  was  treated  as  a  compressible  fluid  in  accordance  with  the 
Ideal  Gas  Law  and  was  subjected  to  adiabatic  constraints  during  bending.  The  fabric  was 
represented  with  membrane  elements.  Several  constitutive  cases,  including  linear  elasticity  and 
hyperelasticity,  were  studied.  The  bending  behavior  for  each  constitutive  case  is  presented, 
followed  by  a  discussion  of  the  use  and  limitations  of  the  cases. 


DESCRIPTION  OF  PLAIN-WOVEN  FABRIC 


The  fabric  was  a  dense,  high-quality,  uncoated,  plain-woven  polyester  as  shown  in  figure  1 
and  is  commercially  used  in  the  construction  of  fire  hoses.  The  warp  yams  were  aligned  in  the 
longitudinal  direction  of  the  fire  hose  and  the  weft  yams,  orthogonal  to  the  warp  yams,  were 
aligned  in  the  hoop  direction.  Specific  details  at  the  yam  and  fabric  levels  were  measured,  using 
reverse  engineering,  and  are  reported  in  table  1 . 
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Figure  1.  Image  of  Plain-Woven  Polyester  Fabric 


Table  1.  Details  of  the  Plain-Woven  Polyester  at  the  Yarn  and  Fabric  Levels 


Characteristic 

Description 

Warp  Yams 

Polyester 

5  bundles 

2  twists/in. 

Cross-section  area  =  0.00105  in.2 

Weft  Yams 

Polyester 

1  bundle 

2  twists/in. 

Cross-section  area  =  0.00153  in.2 

Fabric  Densities 

Warp  yarns  per  unit  circumference  =  15  yarns/in. 

Weft  yams  per  unit  length  of  air  beam  =  1 1  yams/in. 
Yam  density  ratio  (YDR)  =  0.73 

Crimp  Measurements 

cwarp=  15.0% 

Cw/,  =  3.1% 

<T=  Cwarp/Cwef,  =  0.2 1 

The  crimp  ratio  C  is  defined  as  the  waviness  of  the  yams  as  obtained  by  equation  (1): 

^ yarn  ^ fabric  (  | 

L fabric 

where  Lfabnc  is  the  measured  length  of  a  yam  in  its  fabric  state  and  Lyarn  is  the  length  of  that  same 
yam  after  it  has  been  extracted  from  the  fabric  and  straightened  out.  It  was  particularly  noted 
that  the  fabric  was  engineered  with  different  crimp  amounts  between  the  two  yam  families.  The 
crimp  content  of  the  warp  yams  was  nearly  five  times  greater  than  that  of  the  weft  yams.  The 
yam  density  ratio  YDR  was  defined  as  the  number  of  weft  (hoop)  yams  per  unit  beam  length 
divided  by  the  number  of  warp  (longitudinal)  yams  per  unit  circumference. 
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FABRIC  STRESS-STRAIN  BEHAVIOR— EXPERIMENT 


For  pressurized  cylindrical  beams  of  radius  r,  it  can  readily  be  shown  from  equilibrium  that 
the  stress  ratio  S  (defined  as  the  ratio  of  hoop  stress  per  unit  beam  length  to  the  longitudinal 
stress  per  unit  circumference)  is  2: 1,  as  shown  in  figure  2.  Therefore,  to  properly  characterize 
the  stress-strain  behavior  of  the  fabric  about  a  given  pressurized  state,  fabric  samples  must  be 
tested  under  biaxial  loading  conditions  that  enforce  S  =  2.  This  ratio  holds  only  for  the  inflated 
state  and  is  not  applicable  to  pressurized  cylinders  subjected  to  external  loads. 

The  authors  did  not  have  access  to  a  combined  biaxial  and  shear  test  fixture  capable  of  S  =  2 
loading  at  the  time  of  this  research;  the  equi-biaxial  extension  test  fixture3  shown  in  figure  3  was 
used.  For  the  woven  air  beam  models  evaluated  in  this  report,  the  longitudinal  axis  was 
coincident  with  the  warp  direction  and  the  hoop  axis  was  coincident  with  the  weft  direction 
(refer  to  figure  1). 


Weft  tensile  stress  per  unit  length  of  cylinder  -  Pr 
Warp  tensile  stress  per  unit  length  circumference  =  Pr  2 


Figure  2.  Stress  State  in  Cylinders  Subjected  to  Internal  Pressure 


Biaxial  Tension  Mode  In-Plane  Shear  Mode 

Figure  3.  Combined  Multi-Axial  and  Shear  Test  Fixture  (U.S.  Patent  No.  6,860,156) 
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A  swatch  of  the  polyester  fabric,  of  the  dimensions  given  in  figure  4,  was  subjected  to  a 
series  of  uniaxial  tension,  biaxial  tension,  and  combined  biaxial  tension  and  shear  tests,  as 
illustrated  in  figure  5. 


Figure  4.  Dimensions  of  Fabric  Swatch  Specimen 


Figure  5.  Plain-Woven  Polyester  Fabric  Specimen  Subjected  to  Multi-Axial  and  Shear  Tests 

Five  experimental  stress-strain  curves  were  obtained.  These  curves  corresponded  to 
(1)  uniaxial  tension  along  the  warp  direction,  (2)  uniaxial  tension  along  the  weft  direction, 

(3)  biaxial  tension  along  the  warp  direction,  (4)  biaxial  tension  along  the  weft  direction,  and 
(5)  planar  shear. 

The  biaxial  curves — items  3  and  4  listed  above — were  obtained  from  a  single  test  in  which 
the  fixture  operated  in  an  equi-biaxial  extension  rather  than  in  the  equi-biaxial  force  mode.  The 
results  of  the  uniaxial  and  biaxial  tests — items  1  through  4  listed  above — are  plotted  in  figure  6. 
Because  the  maximum  fixture  load  was  limited,  the  fabric  was  not  tested  to  failure. 
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Figure  6.  Uniaxial  and  Biaxial  Stress-Strain  Curves  for  Plain-Woven  Polyester  Fabric 


The  variation  of  S  resulting  from  the  equi-biaxial  extension  mode  of  testing  is  plotted  in 
figure  7.  The  result  is  neither  controlled  nor  constant  because  the  elastic  moduli  of  the  fabric 
were  highly  nonlinear  with  respect  to  extension  (and  strain).  At  the  maximum  equi-biaxial  strain 
e  applied  during  testing,  S  is  3.82,  nearly  twice  the  required  value  for  properly  characterizing  the 
state  of  stresses  in  a  pressurized  cylinder.  However,  the  effects  of  crimp  interchange  during  the 
equi-biaxial  extension  test  clearly  resulted  in  pronounced  differences  between  the  uniaxial  and 
biaxial  curves  for  a  given  yam  family.  Figure  8  shows  the  strain-dependent  elastic  modulus, 
now  referred  to  as  the  tangent  modulus  Elan,  along  the  warp  direction  as  a  function  of  warp  stress 

&warp‘ 


Equi-Biaxial  Strain  (in./in.) 


Figure  7.  Stress  Ratio  S  vs  Equi-Biaxial  Strain  e 
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Figure  8.  E,an  vs  owarp  Measured  in  the  Equi-Biaxial  Extension  Test 


The  planar  shear  test  was  conducted  in  a  constant  equi-biaxial  extension  mode,  and  owarp 
corresponded  to  that  of  a  10-psi  inflation  pressure — a  stress  of  637  psi,  as  shown  in  figure  8.  The 
specimen  was  loaded  to  a  maximum  shearing  angle  of  10°,  and  the  resulting  curve  for  shear 
stress  r  as  a  function  of  shear  strain  y  is  plotted  in  figure  9. 

The  shear  modulus  G,  along  path  ab  in  the  direction  of  increasing  y,  was  derived  by 
differentiating  a  polynomial  curve  fit  of  figure  9.  The  resulting  shear  modulus  is  plotted  in 
figure  10.  Figure  1 1  illustrates  the  three  distinct  regions  of  shear  stiffness,  revealed  by  figure  9, 
that  resulted  from  kinematic  interactions  and  scissoring  deformations  between  yam  families. 


Figure  9.  t  vs  y  Plotted  with  P  at  10  psi  (owarp  =  637  psi) 
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Figure  10.  G  vs  y  with  P  at  10  psi  (awarp  =  637  psi) 


Region  I 

(Onset  of  Slip) 


Region  II 

(Scissoring) 


Region  III 

(Jamming) 


i  ii 
i  i  i 


Figure  11.  Kinematics  of  Shear  Deformations  in  Uncoated  Plain-Woven  Fabrics 


Referring  to  region  I  in  figure  9,  G  was  initially  governed  by  compaction  forces  between  the 
yams  at  the  crossover  points  because  of  the  biaxial  tensile  stresses.  As  y  increased,  yam  slip 
initiated,  and  at  y  equal  to  approximately  0.03  in./in.,  G  declined  to  practically  zero.  At  this 
strain,  rotations  between  yam  families  were  met  with  no  appreciable  resistance.  Upon  further 
increase  of  y  into  region  II,  G  increased  as  the  gaps  between  yam  families  diminished,  and  the 
onset  of  shear  jamming  occurred. 

Region  III  is  known  as  the  shear  jamming  state.  The  onset  of  shear  jamming,  which  can  be 
determined  through  geometric  models,5  is  related  to  the  maximum  number  of  weft  yams  that  can 
be  woven  into  the  fabric  for  a  given  warp  yarn  size  and  spacing.  The  shearing  behavior  of  this 
fabric  was  consistent  with  the  numerical  model  and  test  results  of  the  plain-woven  fabric 
investigated  in  Cavallaro  et  al.2 
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AIR  BEAM  FINITE  ELEMENT  MODELS 


The  air  beams  modeled  in  the  effort  described  in  this  report  were  98  inches  long,  had  a 
nominal  4.0-inch  diameter,  and  were  subjected  to  four-point  bending  as  shown  in  figure  12. 
Although  a  thin  urethane  bladder  is  typically  used  to  line  the  inner  surface  of  air  beams  to 
prevent  air  leakage  through  the  fabric,  it  is  considered  nonstructural.  Therefore,  such  a  bladder 
was  not  included  in  the  models  described  here.  The  models  were  developed  with  Altair’s 
HyperWorks6  pre-  and  post-processor  and  the  ABAQUS/Explicit  solver.4  This  solver  was  used 
to  capture  the  pressure-volume  behavior  of  the  internal  air,  transverse  shear  deformations  from 
bending,  geometric  nonlinearities  due  to  large  deformations,  material  nonlinearities  (for  the 
hyperelastic  material  models),  localized  wrinkling,  and  surface-to-surface  contact  kinematics. 

The  fabric  was  idealized  as  being  incapable  of  developing  bending  strain  energy  and  was 
discretized  using  membrane  elements.  However,  strain  energies  associated  with  extension  and 
shearing  deformations  were  captured.  The  computation  of  the  membrane  element  thickness  was 
based  on  an  analogous  homogeneous  cylinder  with  an  outer  diameter  of  4.00  inches  and  a  cross- 
sectional  area  equivalent  to  the  sum  of  the  warp  (axial)  yam  areas.  The  resulting  thickness  was 
0.016  inch.  Rigid  saddles  that  were  4.00  inches  long  and  4.27  inches  in  inner  diameter,  as  shown 
in  figure  12,  were  used  as  the  load  points  (upper  saddles)  and  support  points  (lower  saddles). 

Full  surface-to-surface  contact  definitions  were  included  between  the  fabric  and  saddle  elements. 


'mid 


Figure  12.  Air  Beam  Finite  Element  Model  Subjected  to  Four-Point  Bending 


The  models  were  loaded  by  enforcing  a  lateral  displacement  Sioad_pi  of  6.00  inches  (1.5  times 
the  diameter)  at  both  load  points.  The  support  points  were  restricted  from  translating  in  all 
directions.  All  the  saddles  were  allowed  to  rotate  so  that  the  lateral  reaction  force  vectors 
remained  normal  to  the  air  beam  at  all  times  (i.e.,  follower  mode).  A  zero  coefficient  of  friction 
was  used  at  all  contact  surfaces  so  that  relative  slip  between  the  air  beam  and  saddles  occurred 
without  restriction.  The  distance  between  support  point  centerlines  was  72.0  inches,  and  the 
distance  between  loading  point  centerlines  was  37.0  inches. 
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To  model  the  fluid-structure  interactions  of  the  internal  air  and  surrounding  fabric,  a 
pressurized  cavity  was  defined  along  the  inside  surface  of  the  fabric  material.  The  cavity  and  its 
enclosing  surface  were  used  to  apply  the  internal  pressure  P  directly  to  the  membrane  (fabric) 
elements  and  to  define  the  volume  of  air  contained  by  the  cavity  within  the  air  beam.  The 
internal  air  was  modeled  as  a  compressible  (pneumatic)  fluid  that  satisfied  the  ideal  gas  equation 
of  state  (EOS)  as  described  in  appendix  A.  This  EOS  assumed  that  compressibility  of  the  air 
occurred  adiabatically;  that  is,  no  heat  transfer  was  permitted  across  the  boundaries  between  the 
cavity  and  the  fabric. 

The  external  work  Wex,  done  on  the  air  beam  during  a  quasi-static  four-point  bend  test  was 
computed  as  the  area  under  the  total  reaction  force  Freaa  as  a  function  of  the  enforced  load  point 
displacement  8joaci_p,  curve.  During  the  bending  step,  W^,  was  equal  to  the  change  in  the  internal 
energy  of  the  air  beam  AFin,.  The  change  in  internal  energy  consisted  of  the  sum  of  the  changes 
in  the  fabric  strain  energy  A Estram,  the  kinetic  energy  of  the  total  system  mass  A Emetic,  the  work 
done  by  compressing  the  air  A\PdV,  and  the  viscous  dissipation  energy  AEdis  due  to  damping. 
(Note  that  A  refers  to  the  change  between  the  bending  and  inflation  end  states.)  An  energy 
balance  was  performed  as  shown  in  equation  (2): 

^  =  K^d^,=A£,„,,  (2) 

where 

A E  —  A E .  +  AF  +  A\PdV + AE,  . 

Transverse  shear  deformations7  arising  from  the  shearing  strain  y,  as  shown  in  the  beam¬ 
bending  example  in  figure  13,  will  cause  a  reduction  in  air  volume  and  lead  to  additional 
PdV-work. 


Undeformed  State 


Neutral 

Plane 


Figure  13.  Transverse  Shear  Deformations  During  Bending 
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Although  the  ABAQUS/Explicit  solver  was  formulated  for  dynamic  problems,  it  can  be  used 
for  static  applications  as  well.  The  analyst  must  ensure  that  no  modes  of  the  structure  are  excited 
for  static  problems.  One  method  is  to  ensure  that  loads  are  applied  slowly  and  that  A Ekinelic  is 
small.  However,  if  wrinkling  or  buckling  events  were  to  occur,  spikes  in  the  A Ekmeiic  time  history 
curve  would  develop  and  allow  for  easy  detection  of  the  event  in  time.  For  the  effort  described 
here,  A Ekmeiic  and  A Ejis  were  typically  observed  to  make  negligible  contributions  to  equation  (2). 

To  evaluate  the  bending  performance  of  the  air  beam,  a  two-step  loading  process  was 
necessary.  First,  the  cavity  was  pressurized  to  the  required  inflation  level  during  a  2-second  time 
interval.  Second,  dioad_pi =  6.0  inches  was  applied  at  the  load  points  during  a  5-second  time 
interval. 

The  following  section  describes  the  implementation  of  material-level  test  results  in  the  air 
beam  bending  models.  Specifically,  the  warp  direction  stress-strain  curve  taken  from  the 
equi-biaxial  extension  test  was  used  to  formulate  the  linearly  elastic  and  hyperelastic  constitutive 
behaviors  assigned  to  the  fabric  membrane  elements. 


ANALYSIS  AND  RESULTS 


CASE  1— LINEAR  ELASTICITY,  CONSTANT  MODULUS  FOR  ALL  PRESSURES 

The  first  case  idealized  the  fabric  material  as  isotropic  and  linearly  elastic.  The  elastic  and 
shear  moduli  were,  therefore,  invariant  with  inflation  pressure.  Here,  the  shear  modulus  G  and 
elastic  modulus  E  were  related  through  a  generalized  form  of  Hooke’s  Law  as 


G  = 


E 

2(1  +  y)’ 


(3) 


where  vis  Poisson’s  ratio.  The  value  of  E  was  taken  as  the  largest  tangent  modulus  of  the  warp 
direction  biaxial  stress-strain  curve  E,a„  (0. 1  Mpsi),  and  v  was  set  to  0.3  to  allow  for  material 
compressibility. 

The  time  histories  of  pressure,  internal  volume,  load-point  reaction  forces,  and  energy  terms 
necessary  to  conduct  the  energy  balance  of  equation  (2)  were  tracked  during  the  solution.  The 
change  in  pressure  A P  was  plotted  as  a  function  of  the  corresponding  change  in  volume  AV 
during  bending.  The  area  under  this  curve  was  integrated  to  obtain  PdV-work. 

For  this  linearly  elastic  case  and  the  inflation  pressures  of  interest,  it  was  observed  that  A P  and 
AV  were  very  small  during  bending.  The  results  indicated  that  A Eslrain  decreased  with  increasing 
P,  while  PdV-work  increased  with  increasing  P,  as  shown  in  figure  14  for  the  specified  5ioadj>t-  At 
pressures  below  30  psi,  PdV-work  contributed  less  than  10%  of  W, and  A Es,rai„,  contributing 
more  than  90%,  was  the  dominant  term  in  the  energy  balance.  However,  for  inflation  pressures  of 
90  psi  and  greater,  PdV-work  dominated  the  energy  balance,  almost  to  the  exclusion  of  AEslrain. 
The  effects  of  A Es,rai„  on  bending  behavior  were  negligible  at  these  pressures. 
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Figure  14.  Wexl,  PdV-work,  and  \Eslrain  vs  P  for  a  Linearly  Elastic  Fabric — Case  1 


Figure  15  shows  the  curves  for  the  total  reaction  force  Freac,  versus  8mui  that  were  plotted  for 
the  pressures  of  interest.  For  the  enforced  5i0ad_pi  of  6.0  inches,  the  occurrence  of  wrinkling  was 
predicted  for  pressures  of  10  to  40  psi  only.  The  pressures  shown  in  figure  15  for  the  onset  of 
wrinkling  were  consistent  with  the  theoretical  values  obtained  from  equation  (4), 


M 


wrinkle 


Ptt  r  3 
2 


(4) 


which  is  easily  derived  from  a  simple  stress  balance  between  the  inflated  and  bending  stress 
states  in  woven  air  beams. 


Figure  15.  Freact  vs  6m  -,d  Curves  for  Linear  Elastic  Material  Model  with  P  at  10  to  60  psi 

— Case  1 
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In  comparison  with  previous  air  beam  bending  experiments,  the  curves  in  figure  15  exhibit  a 
moderately  weak  dependence  on  P.  This  effect  is  attributed  to  the  linearly  elastic  constitutive 
relationship  of  equation  (3),  whereby  G,  computed  as  38,462  psi,  is  overly  stiff  in  comparison 
with  the  experimental  results  from  figure  10  that  reveal  a  peak  G  of  4,600  psi,  corresponding  to 
P  at  10  psi.  As  a  result,  no  appreciable  transverse  shearing  deformations  develop  in  these 
models,  and  AV,  A P  and  PdV-work  are  inconsequential  to  the  bending  behavior.  These 
parameters  are  summarized  for  the  bending  step  in  table  2. 


Table  2.  Summary  of  Pressure,  Volume,  and  Energy  Changes  During  Bending  Step — Case  1 


p 

(psi) 

A P 
(psi) 

AV 

(in.3) 

AlPdV 

(in.-Ib) 

AE strain 

(in.-Ib) 

10 

0.017 

-0.88 

8.75 

169.64 

20 

0.024 

-0.89 

18.20 

171.03 

30 

0.019 

-0.56 

16.88 

179.65 

40 

0.025 

-0.63 

25.07 

180.65 

50 

0.036 

-0.79 

39.25 

176.69 

60 

0.055 

-1.07 

64.23 

164.81 

CASE  2— LINEAR  ELASTICITY,  USE  OF  INSTANTANEOUS  TANGENT  MODULI 

This  case  assumed  that  the  fabric  material  behaved  as  linearly  elastic  in  accordance  with 
equation  (3),  but  E,a„  used  at  each  pressure  corresponded  to  the  instantaneous  tangent  modulus 
E ',an  at  the  state  of  stress  because  of  inflation  at  that  pressure.  In  essence,  case  2  allowed  E,an  and 
G  to  change  with  respect  to  the  specific  inflation  pressure  only.  From  the  warp  direction  biaxial 
stress-strain  curve,  the  resulting  values  of  owarp  and  E',an  are  computed  for  inflation  pressures  of 
10,  15,  and  20  psi  (see  table  3).  A  28%  increase  in  E',an  is  observed  by  increasing  the  inflation 
pressure  from  10  psi  to  20  psi. 


Table  3.  owarp  and  E',an  as  Functions  of  P 


P 

(Psi) 

G  warp 

(psi) 

Eft 

c  tan 

(psi) 

10 

635 

58,765 

15 

952 

69,924 

20 

1,270 

75,414 

Case  2  was  demonstrated  for  inflation  pressures  only  up  to  20  psi  because  load  limitations  of 
the  biaxial  test  fixture  prevented  characterization  of  the  biaxial  warp  stress-strain  curve  beyond 
1,500  psi.  Hence,  E',an  could  not  be  determined  for  stresses  beyond  this  limit.  In  addition,  the 
results  from  this  case  were  valid  for  only  those  bending  stresses  that  could  be  treated  as 
perturbations  from  the  inflated  stress  states  where  the  tensile  and  compressive  bending  behaviors 
did  not  deviate  from  E',a„  of  the  inflated  stress  state. 
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The  results  for  the  work  and  energy  terms  were  plotted  as  functions  of  inflation  pressure  in 
figure  16.  The  values  from  case  2  were  similar  to  those  from  case  1.  Curves  of  Freac,  versus  6miet 
were  plotted  in  figure  17  for  P  at  10,  15,  and  20  psi.  For  case  2,  wrinkling  was  predicted  for 
each  pressure  as  shown  in  figure  17.  Changes  in  the  pressure,  volume,  and  energy  terms  were 
summarized  in  table  4. 


Figure  16.  We__xt,  A£s,raim  and  PdV -work  vs  P for  a  Linearly  Elastic  Fabric  When  E',an  Used 

with  P  at  10, 15,  and  20  psi — Case  2 


Figure  1 7.  Freact  vs  SmM  for  a  Linearly  Elastic  Fabric  When  E  Used  with  P  at  10, 15, 

and  20  psi — Case  2 
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Table  4.  Summary  of  Pressure,  Volume,  and  Energy  Changes  During  Bending  Step — Case  2 


p 

(psi) 

A P 
(psi) 

W 

(in.3) 

aJ/w 

(in.-lb) 

A E sirain 

(in.-lb) 

10 

0.017 

-0.86 

8.65 

105.29 

15 

0.017 

-0.73 

10.96 

127.90 

20 

0.016 

-0.62 

12.34 

146.25 

While  further  characterization  of  the  fabric  constitutive  properties  (notably,  the  shear 
modulus)  at  pressures  beyond  20  psi  is  necessary,  the  influences  of  pressure  in  case  2 
qualitatively  followed  those  observed  in  previous  bending  experiments.1,2  As  a  result  of  the 
changes  in  both  Estrmn  and  PdV-work  with  respect  to  P,  plain-woven  air  beams — in  particular — 
must  consider  both  terms.  Plain-woven  air  beams  typically  operate  at  low-pressure  levels  (less 
than  those  for  triaxial-woven  or  braided  air  beams)  because  of  safety  concerns  and  are,  therefore, 
susceptible  to  greater  transverse  shear  deformations. 


CASE  3— HYPERELASTICITY 

This  case  considered  both  geometric  and  material  nonlinearities.  Here,  the  full  nonlinear 
constitutive  behavior  of  the  plain-woven  polyester  fabric  was  idealized  as  a  hyperelastic 
material.  The  notable  difference  between  the  hyperelastic  and  materially  linearly  elastic  cases 
was  that  hyperelasticity  enabled  the  stiffness  of  the  fabric  membrane  to  change  with  stresses  due 
to  inflation  pressure  and  bending  loads. 

A  hyperelastic  strain  energy  potential  was  pursued  because  the  uniaxial  and  biaxial  stress- 
strain  curves  obtained  from  swatch-level  tests  exhibited  nonlinear  stiffening  similar  to  that  of 
hyperelastic  materials.  Several  strain  energy  potentials  were  evaluated  using  ABAQUS/Explicit 
to  determine  the  best  representation  of  the  warp  direction  biaxial  stress-strain  curve.  This  stress- 
strain  curve  specifically  dominated  the  material  behavior  of  the  air  beam  during  bending.  The 
reduced  polynomial  ( N  =  3)  strain  energy  potential,4  shown  in  appendix  B,  provided  the  best  fit 
for  this  stress-strain  curve  and  was  stable  over  the  entire  strain  range  (i.e.,  5%).  The  initial  shear 
modulus  po  was  507.2  psi.  The  results  for  the  work  and  energy  terms  are  plotted  as  functions  of 
inflation  pressure  in  figure  18. 

Curves  of  Freaci  versus  6mi(i  are  plotted  in  figure  19,  which  clearly  shows  dependence  of  the 
bending  behavior  on  P.  However,  wrinkling  is  not  predicted  at  these  pressures  and  8ioad_p, 
because  the  moments  generated  are  less  than  required  by  equation  (4).  The  changes  in  pressure, 
volume,  and  energy  terms  during  the  bending  step  are  listed  in  table  5  as  a  function  of  P. 
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Figure  18.  AE strain,  and  PdV-work  vs  P  for  a  Hyperelastic  Fabric:  Reduced  Polynomial 

(TV  =  3)  Strain  Energy  Potential  Applied  with  P  at  10, 15,  and  20  psi — Case  3 


Figure  19.  Freac,  vs  Smidfor  Hyperelastic  Fabric  Model  with  P  at  10, 15,  and  20  psi — Case  3 

Table  5.  Summary  of  Pressure,  Volume,  and  Energy  Changes  During  Bending  Step  for  the 

Hyperelastic  Fabric  Model — Case  3 


p 

(psi) 

\P 

(psi) 

AV 

(in.3) 

Aj PdV 
(in.-lb) 

AE strain 

(in.-lb) 

10 

0.052 

-2.96 

29.74 

16.81 

15 

0.063 

-3.13 

46.94 

17.54 

20 

0.073 

-3.22 

64.66 

20.10 

Time  history  curves  of  owarp  on  the  outermost  compressive  surface  at  the  mid-span  are  plotted 
in  figure  20  for  the  inflation  and  bending  steps  for  P  at  10,  15,  and  20  psi. 
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Figure  20.  Time  History  Plot  of  owarpfor  Inflation  and  Bending  Steps 
with  P  at  10, 15,  and  20  psi — Case  3 


Plots  of  owarp  and  owef,  versus  warp  and  weft  strain,  respectively,  are  shown  in  figure  2 1  for 
the  inflation  and  bending  steps  at  the  mid-span,  outermost  compressive  surface.  Stresses  due  to 
inflation  are  designated  with  dashed  lines,  and  stresses  due  to  bending  are  designated  with  solid 
lines.  Note  that  S  during  inflation  is  approximately  2: 1,  and  only  minimal  changes  in  awef,  occur 
during  bending,  particularly  for  P  at  1 0  and  15  psi. 

The  treatment  of  the  fabric  constitutive  behavior  as  a  hyperelastic  material  is  limited  to  those 
fabrics  at  stress  levels  in  which  the  actual  fabric  shear  modulus  can  be  sufficiently  represented  by 
the  strain  energy  potential. 


Strain  (in./in.)  Strain  (in./in.)  Strain  (in./in.) 

Figure  21.  Plots  of  owarp  and  oweft  v.v  Strain  for  P  at  10, 15,  and  20  psi — Case  3 
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DISCUSSION 


The  bending  performance  of  plain-woven  fabric  air  beams  was  investigated  through  a 
combined  material-level  test  and  finite  element  analysis  approach.  Using  ABAQUS/Explicit,  the 
fluid-structure  interactions  resulting  from  pressure-volume  changes  during  inflation  and 
four-point  bending  events  were  determined.  The  warp  tensile  stress-strain  curve  of  the  fabric 
was  established  from  an  equi-biaxial  extension  test  and  was  used  to  formulate  both  linear  elastic 
and  hyperelastic  constitutive  cases  for  the  air  beam  bending  model. 


The  material-level  tests  confirmed  that  the  plain-woven  fabric  did  not  behave  as  a  continuum 
but  rather  as  a  discrete  assemblage  of  fibers  and  yams.  For  the  range  of  pressures  considered,  E 
and  G  were  independent  of  each  other  and  were  based  on  changes  in  fabric  architecture  (crimp 
interchange,  yam  slip,  rotation,  etc.),  YDR,  and  stress  ratio  S.  Although  E  obtained  from  the 
biaxial  warp  stress-strain  curve  was  highly  nonlinear  with  P,  it  was  monotonic  with  stress  (and 
pressure)  beyond  0.5  psi  and  was  easily  represented  using  a  polynomial  curve  fit  for  input  to  the 
air  beam  models.  However,  G  was  not  monotonic  in  y.  Three  distinct  regions  developed  in  the 
plot  of  r  versus  y  and  were  dominated  by  yam  rotations  and  shear  jamming.  Referring  to  figure 
9,  regions  I  and  II  were  highly  dependent  on  y;  and  region  II  was  expected  to  produce  the  largest 
A  Fin  the  plain-woven  air  beam  models  because  of  transverse  shearing  deformations  from 
bending.  During  bending,  this  shear-induced  AF  was  a  source  of  PdV-work. 

The  material-level  biaxial  test  was  conducted  in  an  equi-biaxial  extension  mode.  As  a 
consequence,  S  was  not  controlled,  but  it  is  hypothesized  that,  for  an  uncoated,  plain-woven 
fabric  subjected  to  equi-biaxial  extension,  the  value  of  S  for  which  convergence  in  e  occurs, 
denoted  S°,  appears  to  be  given  by  equation  (5): 


sc=f\ 


\ 


,C'JDR 

J 


(5) 


To  properly  simulate  the  inflated  air  beam  stress  state  in  a  swatch  specimen,  a  combined  biaxial 
and  shear  test  fixture  with  the  capacity  to  regulate  a  2: 1  biaxial  tensile  stress  ratio  is  needed. 
Development  efforts  are  underway  to  provide  this  experimental  capability. 

Comparisons  were  made  between  the  fabric  strain  energy  Estrai„  and  PdV-work  done  on  the 
air  beam  during  bending.  It  was  shown  in  linear  elastic  case  1  that  Estram  dominated  the  bending 
behavior  when  P  <  30  psi,  so  that  AF  and  PdV-work  were  minimal.  However,  when  P  >  90  psi 
and  AF  was  large,  PdV-work  dictated  bending  behavior,  and  the  material  constitutive  effects 
{E strain)  on  bending  behavior  were  minimal.  Results  from  linear  elastic  case  2  obtained  for 
pressures  up  to  20  psi  also  showed  that  Es,rain  dominated  the  bending  behavior.  The  hyperelastic 
case  3  results  exhibited  larger  AF  and  PdV-work  for  the  pressures  of  10,  15,  and  20  psi. 
However,  Eslraj„  contributed  approximately  25%  to  the  energy  balance.  In  general,  both  terms — 
Estrain  and  PdV-work — must  be  determined  for  air  beams  subject  to  bending  loads  because  y  will 
increase  with  decreasing  pressure  and  lead  to  increasing  values  of  AF. 
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The  elastic  and  hyperelastic  constitutive  theories  are  of  limited  use  if  shearing  deformations 
must  be  considered.  These  theories  couple  the  extension  and  shear  behaviors,  which,  on  the 
contrary,  were  shown  uncoupled  through  material-level  tests.  Consequently,  the  use  of  these 
theories  is  unlikely  to  capture  the  critical  shear  response  (regions  I  and  II  in  figure  9)  and  may 
lead  to  erroneous  bending  predictions.  However,  these  theories  may  provide  acceptable 
constitutive  predictions  for  use  in  models  of  coated  woven  fabric  air  beams.  Aside  from  the 
obvious  environmental  protections  afforded  by  coatings,  coatings  provide  additional  elastic  and 
shear  stiffness  by  minimizing  yam  slip  and  rotation.  The  added  shear  stiffness  reduces  shearing 
deformations,  especially  those  developing  in  regions  I  and  II.  As  a  result,  A  Land  PdV-work  will 
be  less  in  comparison  to  the  same  terms  for  uncoated  fabric  air  beams.  In  the  limiting  case,  the 
constitutive  behavior  of  coated  woven  fabrics  approaches  that  of  matrix-reinforced  fibrous 
composites. 

If  the  fabric  constitutive  properties  are  nonlinear,  the  instantaneous  tangent  modulus,  E\an, 
must  be  used  at  the  corresponding  stress  state  arising  from  the  specific  inflation  pressure. 
However,  if  E',a„  and  G  change  considerably  with  respect  to  pressure  and  bending  loads,  these 
moduli  must  be  numerically  updated  in  the  models  to  properly  assess  bending  behavior  beyond 
simple  perturbation  states.  Hence,  the  need  exists  to  measure  AV  during  air  beam  bending  tests 
so  that  sufficient  model  validation  can  be  performed. 


CONCLUSIONS 


It  is  concluded  that,  for  a  specific  air  beam  deflection,  the  volume  change  will  be  larger  at 
higher  inflation  pressures  than  for  lower  inflation  pressures  because  the  beam  is  stiffer  at  higher 
pressures  and  will  require  more  external  force  to  achieve  the  deflection.  The  higher  external 
force  creates  larger  shear  forces  in  the  beam,  and,  thus,  a  greater  shear  distortion  occurs. 
Therefore,  volume  changes  are  expectedly  greater  with  increasing  inflation  pressures  for  a 
prescribed  deflection. 

Air  compressibility  ( PdV-work )  introduces  a  nonlinear  stiffening  effect  in  the  bending 
behavior  of  inflatable  fabric  structures.  As  the  air  volume  decreases  due  to  deformations  from 
external  loads  (such  as  transverse  shear,  wrinkling,  and  section  collapse),  the  air  pressure  will 
increase.  The  air  behaves  as  a  nonlinear  spring,  and  its  impact  on  bending  behavior  is  directly 
related  to  AV.  If  air  compressibility  is  appreciable,  a  gas  law  must  be  incorporated  in  the 
bending  analysis. 
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APPENDIX  A 

IDEAL  GAS  EQUATION  OF  STATE  (EOS) 


The  ideal  gas  equation  of  state  (EOS)*  is  given  as 

P  +  Pa=Pr{6-0z ), 

where 

P  is  the  internal  pressure, 

Pa  is  the  ambient  pressure  (14.7  psi), 
p  is  the  density  of  air  (4.4274e-005  lb/in.3), 

R  is  the  gas  constant, 

0  is  the  current  temperature  (2 1  °C),  and 
6Z  is  the  absolute  zero  temperature  (-273.15°C). 
The  ideal  gas  constant  R  is  given  by 
R  =  R/MW, 


where 


R  is  the  universal  gas  constant 


2.841  xlO4 


V 


lb  in.2  A 
mole  s2°K 


,  and 


Mw  is  the  molecular  weight  (28.97  moles). 


(A- 1 ) 


(A-2) 


*“ABAQUS/Explicit,”  Version  6.4,  ABAQUS  Inc.,  Pawtucket,  RJ,  2003. 
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APPENDIX  B 

REDUCED  POLYNOMIAL  (N=  3)  STRAIN  ENERGY  POTENTIAL 


The  reduced  polynomial  (N=  3)  strain  energy  potential*  is  given  as 

i= I  /= 1  ui 


(B-l) 


where 

U  =  strain  energy  per  unit  volume, 

N  =  material  parameter, 

C, o  =  temperature  dependent  material  parameter, 

D,  =  temperature  dependent  material  parameter, 

7,  =  first  deviatoric  strain  invariant,  and 

Jel  =  elastic  volume  ratio. 

The  initial  shear  modulus  /.i0  is  given  by 

A,  =2  (CI0  +  C0I).  (B-2) 

To  define  the  limiting  case,  the  fabric  membrane  elements  were  assumed  to  be  fully 
incompressible  by  setting  D\,  Di,  and  D3  to  zero.  The  initial  shear  modulus,  /uo,  was  507.2  psi. 
The  constants  Cn,  C10,  C20,  and  C30  are  given  below: 

C/i  =  C/2  =  0  psi, 

C 10  =  253.62  psi 

C20  =  155,072.53  psi,  and 

C30  =  1,914,968.40  psi. 


*“ABAQUS/Explicit,”  Version  6.4,  ABAQUS  Inc.,  Pawtucket,  RJ,  2003. 
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